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^^ ' In this talk we discuss some details of the analysis of neutrino data and our present 

C 3 ' understanding of neutrino masses and mixing. This talk is based on Refs. [1, 2, 3]. 



1. Analysis of Solar and KamLAND 



~~f • In Fig. 1 we sliow tlie results from our latest analysis of KamLAND [4] 

reactor 17^ disappearance data, solar Vg data [5, 6] and the combined analysis 



under the hypothesis of CPT symmetry. The main new ingredient in this 
analysis with respect to the previous ones is the inclusion of the first results 
from the SNO salt phase (SNOII) data [6]. We have also taken into account 
the new gallium measurement which leads to the new average value 69.3 ± 
4.0. The main changes as compared to the pre- SNOII analysis are: 

- in the analysis of solar data, only LMA is allowed at more than 3ct; 

- maximal mixing is rejected by the solar analysis at more than Scr; 

- the combined analysis allows only the lowest LMA region at 99.4% CL; 

- the new best-fit point is: 

Am^ = 7.1 X 10"^ eV^ tan^ 6l = 0.41 , '^ ^0.88. (1) 

't'8B,BP04 

These results are in agreement with those reported in the several state-of- 
the-art analysis of solar and KamLAND data which exist in the literature. 
All these analysis share the same basic characteristics. 

In the analysis of KamLAND we used the following approximations: 

- the antineutrino spectrum is parameterized [8] without detailed theoret- 
ical uncertainties; 
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Figure 1. Allowed regions for 2-u oscillations of Pe in KamLAND and CHOOZ (left 
panel), of v^ in the sun (central panel), and for the combination of KamLAND and 
solar data under the hypothesis of CPT conservation (right). The different contours 
correspond to the allowed regions at 90%, 95%, 99% and 3cr CL. 



- the yearly average reactor power is used. 

Presently the uncertainties of the KamLAND results are statistics domi- 
nated so these effects do not make any difference in the extracted allowed 
regions. The minor differences between the several phenomenological anal- 
ysis in the literature are more likely to arise from the use of different sta- 
tistical functions in the analysis of KamLAND data. 

Before moving to atmospheric neutrinos, we wish to point out some 
important features of the analysis of solar data: 

- the SSM [7] provides detailed informations not only on the solar neu- 
trino fluxes themselves, but also on their theoretical uncertainties and 
correlations due to variations of the SSM inputs; 

- the spectral shape experimental uncertainty for ^B spectrum is properly 
taken into account; 

- the energy dependence of the interaction cross sections and their uncer- 
tainties is also included; 

- the interplay between the energy-dependent part of the theoretical and 
systematic uncertainties and the neutrino survival probability (which 
depends on the oscillation parameters) is properly taken into account. 



2. Atmospheric Neutrinos 

In the left panel of Fig. 2 we show the results of our latest analysis of 
the atmospheric neutrino data, which included the full data set of Super- 
Kamiokande phase I (SKI). As discussed in Ref. [9] the new elements in 
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Figure 2. Le/4.' allowed regions from the analysis of atmospheric data using the new 
(full regions labeled "3-dim") and old (empty curves labeled as "1-dim") SKI data and 
atmospheric fluxes. Right: impact of a small change in the correlation between the 
theoretical uncertainties of the low-energy (sub-GeV) and high-energy (multi-GeV, and 
upgoing fi) data samples. The different contours correspond to at 90%, 95%, 99% and 
So- CL. 



the Super-Kamiokande analysis include: 

- use of new three-dimensional fluxes from Honda [fO]; 

- improved interaction cross sections which agree better with the measure- 
ments performed with near detector in K2K [11]; 

- some improvements in the Monte-Carlo which lead to some changes in 
the actual values of the data points. 

We have included these elements in our calculations and we have also im- 
proved our statistical analysis (see Ref. [3] for details). Our results show 
good quantitative agreement with those of the Super-K collaboration. In 
particular we find that after inclusion of the above effects, the allowed re- 
gion is shifted to lower Arn^. The new best-fit point is located at: 



Am^ ^ 2.2 X 10"^ eV^ 



sin^ e = 0.5. 



(2) 



At this point it is important to remark that, unlike for solar neutrinos, 
the energy dependence of the theoretical uncertainties in the atmospheric 
fluxes and in the interaction cross sections are not so well determined in 
terms of a set of model inputs. To address this issue, we have performed the 
same analysis assuming slightly different correlations among the theoretical 
errors for the different data sets, and we have found that the size of the 
final shift in Ani^ of the allowed region depends on these details (see right 
panel of Fig. 2). The reason for this behavior is illustrated in Fig. 3, where 
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Figure 3. The left (right) panels show the allowed regions for the analysis of the different 
sets of atmospheric data using the new (old) SKI analysis and atmospheric fluxes. The 
different contours correspond to the allowed regions at 90%, 95%, 99% and 3cr CL. 



we show the ahowed regions obtained with the new analysis as compared 
to the old one for the different atmospheric data sets. As can be seen in 
the figure, the different sets favor slightly different ranges of Ani^, thus the 
treatment of the energy dependence of the uncertainties becomes relevant 
in outcome of the combined analysis. 

These results lead us to raise here a word of caution. In all present anal- 
ysis of atmospheric data, two main sources of theoretical flux uncertainties 
are included: an energy independent normalization error and a "tilt" er- 
ror which parametrizes the uncertainty in the E~'^ dependence of the flux. 
Some additional uncertainties in the ratios of the samples at different en- 
ergies are also allowed as well as uncertainties in the zenith dependences. 
However, we still lack a well established range of theoretical flux uncer- 
tainties within a given atmospheric flux calculation, in a similar fashion to 
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what it is provided for the solar neutrino fluxes by the SSM. In the absence 
of these, we cannot be sure that we are accounting for the most general 
characterization of the energy dependence of the atmospheric neutrino flux 
uncertainties. 

Given the large amount of data points provided by the Super-K ex- 
periment, this is becoming an important issue in the atmospheric neutrino 
analysis. There is a chance that the atmospheric fluxes may be still too 
"rigid" , even when allowed to change within the presently considered un- 
certainties. As a consequence, we may be over-constraining the oscillation 
parameters. 

3. Three-Neutrino Oscillations 

The minimum joint description of atmospheric [9], K2K [11], solar [5, 6] 
and reactor [4, 12] data requires that all the three known neutrinos take 
part in the oscillation process. The mixing parameters are encoded in the 
3x3 lepton mixing matrix which can be conveniently parametrized in the 
standard form: 

C21 Si2 0\ 
[/ = I C23 S23 I X ( 1 I X ( -Si2 Ci2 (3) 

1/ 

where Cij = cos 9ij and Sij = sin 9ij. Note that the two Majorana phases 
are not included in the expression above since they do not affect neutrino 
oscillations. The angles 9ij can be taken without loss of generality to lie in 
the first quadrant, % S [0,7r/2]. 

There are two possible mass orderings, which we denote as "normal" 
and "inverted" . In the normal scheme rrii < ?7i2 < ms while in the inverted 
one ms < mi < m2- The two orderings are often referred to in terms of 
sign(ATO|i). 

In total the three-neutrino oscillation analysis involves seven parame- 
ters: 2 mass differences, 3 mixing angles, the CP phase and the sign of 
Ato§]^. Generic three-neutrino oscillation effects include: 

- coupled oscillations with two different oscillation lengths; 

- CP violating effects; 

- difference between Normal and Inverted schemes. 

The strength of these effects is controlled by the values of the ratio of mass 
differences a = Ato2i/|A?7i|]^|, by the mixing angle ^13 and by the CP 
phase 5. 
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For solar and atmospheric oscillations, the required mass differences 
satisfy: 

Ami = ^^21 « l^^lil = ^^Itm- (4) 

Under this condition, the joint three-neutrino analysis simplifies and we 
have: 

- for solar and KamLAND neutrinos, the oscillations with the atmospheric 
oscillation length are completely averaged and the survival probability 
takes the form: 

P!^ = Sin^ 013 + cos^ 013P2- (5) 

where in the Sun P^^ is obtained with the modified sun density N^ — » 
cos^^iaA^'g. So the analyses of solar data constrain three of the seven 
parameters: A?7i2i,6'i2 and ^13. The effect of ^13 is to decrease the 
energy dependence of the solar survival probability; 

- for atmospheric and K2K neutrinos, the solar wavelength is too long 
and the corresponding oscillating phase is negligible. As a consequence, 
the atmospheric data analysis restricts Am^i ~ AtoIj, ^23 and ^13, the 
latter being the only parameter common to both solar and atmospheric 
neutrino oscillations and which may potentially allow for some mutual 
influence. The effect of ^13 is to add a ly^ ^ i^e contribution to the 
atmospheric oscillations; 

- at CHOOZ [12] the solar wavelength is unobservable if Arn^ < 8 x 
f 0^** eV and the relevant survival probability oscillates with wavelength 
determined by Am'^i and amplitude determined by ^13. 

In this approximation, the CP phase is unobservable. In principle there 
is a dependence on the Normal versus Inverted orderings due to matter 
effects in the Earth for atmospheric neutrinos. However, this effect is con- 
trolled by the mixing angle ^13, which is constrained to be small by the 
combined analysis of CHOOZ reactor and atmospheric analysis. As a con- 
sequence, this effect is too small to be statistically meaningful in the present 
analysis. 

In Fig. 4 we plot the individual bounds on each of the five parame- 
ters derived from the global analysis. To illustrate the impact of SNOII 
and the new ATM analysis we also show the corresponding bounds when 
either SNOII and the new ATM fluxes are not included in the analysis. 
In each panel the displayed x^ has been marginalized with respect to the 
undisplayed parameters. 
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Figure 4. Global 3u oscillation analysis. Each panel on the left shows the dependence 
of Ax^ on each of the five parameters from the global analysis (full line) compared to 
the bound prior to the inclusion of the new ATM analysis (dashed blue line) and SNO 
data (dotted black line). 



Quantitatively we find the following 3cr CL allowed ranges: 



5.2 < Aml^/IQ-^ eV^ < 9.8, 
1.4 < Aml^/lO^^ eV^ < 3.4, 



0.29 < tan^ 6I12 < 0.64 , 
0.49 < tan^ 6*23 < 2.2 , 

sin^ 613 < 0.054 . 



(6) 



These results can be translated into our present knowledge of the moduli 
of the mixing matrix U: 



\U\ 
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(7) 
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3.1. LSND and Sterile Neutrinos 

Together with the results from the solar and atmospheric neutrino experi- 
ments, we have one more piece of evidence pointing towards the existence of 
neutrino masses and mixing: the LSND experiment, which found evidence 
of F^ -^ Ve neutrino conversion with Aw? > 0.1 eV^. All these data can 
be accommodated into a single neutrino oscillation framework only if there 
are at least three different scales of neutrino mass-squared differences. This 
requires the existence of a fourth light neutrino, which must be sterile in 
order not to affect the invisible Z° decay width, precisely measured at LEP. 

One of the most important issues in the context of four-neutrino sce- 
narios is the neutrino mass spectrum. There are six possible four-neutrino 
schemes which can in principle accommodate the results of solar and at- 
mospheric neutrino experiments as well as the LSND result. They can be 
divided in two classes: (3-1-1) and (2-1-2). In the (3+1) schemes, there is a 
group of three close-by neutrino masses that is separated from the fourth 
one by a gap of the order of 1 eV , which is responsible for the SBL os- 
cillations observed in the LSND experiment. In (2-1-2) schemes, there are 
two pairs of close masses separated by the LSND gap. The main difference 
between these two classes is the following: if a (2+2)-spectrum is realized 
in nature, the transition into the sterile neutrino is a solution of either the 
solar or the atmospheric neutrino problem, or the sterile neutrino takes 
part in both. This is not the case for a (3+l)-spectrum, where the ster- 
ile neutrino could be only slightly mixed with the active ones and mainly 
provide a description of the LSND result. 

The phenomenological situation at present is that none of the four- 
neutrino scenarios are favored by the data. Concerning (2-|~2)-spectra, 
they are ruled out by the existing constraints from the sterile oscillations 
in solar and atmospheric data. As for (3+l)-spectra, they are disfavored 
by the incompatibility between the LSND signal and and the negative re- 
sults found by other short-baseline laboratory experiments. There is also 



February 2, 2008 0:2 WSPC/Trim Size: 9in x 6in for Proceedings 



noon04 




sin 2eL,^o 



0.2 0.4 O.f 



Figure 5. Left: Status of the 3+1 oscillation scenarios. Right: Present status of the 
bounds on the active-sterile admixture from solar and atmospheric neutrino data in 
(2+2)-models. 



a constraint on the possible value of the heavier neutrino mass in this sce- 
nario from their contribution to the energy density in the Universe which is 
presently constrained by cosmic microwave background radiation and large 
scale structure formation data [14]. 

We show in Fig. 5 the latest results of the analysis of neutrino data 
in these scenarios. In the left and central panel we summarize the results 
from Ref. [13] on the (3-1-1) scenarios; we see that after the inclusion of 
the cosmological bound there is only a marginal overlap at 95% CL be- 
tween the allowed LSND region and the excluded region from SBL-I-ATM 
experiments. The right panel illustrates the status of the (2-1-2) scenarios. 
At present, the lower bound on the sterile component from the analysis of 
atmospheric data and the upper bound from the analysis of solar data do 
not overlap at more than Aa. The figure also illustrates the effect of the 
inclusion of the SNOII in this conclusion. 

Alternative explanations to the LSND result include the possibility of 
CPT violation [15], which implies that the masses and mixing angles of neu- 
trinos may be different from those of antineutrinos. We have performed an 
analysis of the existing data from solar, atmospheric, long baseline, reactor 
and short baseline data in the framework of CPT violating oscillations [2] . 
The summary of the results of this analysis is presented in Fig. 6, which 
shows clearly that there is no overlap below the 3cr level between the LSND 
and the all-but-LSND allowed regions. We also note that that the all-but- 
LSND region is restricted to Arn^i = Attilsnd < O-^^ eV^, whereas for 
LSND we always have Arnl^ = Amlg^Y) > 0.02 eV^. 
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